Charge transport through a SET with a mechanically oscillating island 
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We consider a single-electron transistor (SET) whose central island is a nanomechanical oscillator. 
The gate capacitance of the SET depends on the mechanical displacement, thus, the vibrations of 
the island vibrations may strongly influence the current-voltage characteristics, current noise, and 
higher cumulants of the current. Ifarmonic oscillations of the island and oscillations with random 
amplitude (e.g., due to the thermal activation) change the transport characteristics in a different 
way. The noise spectrum has a peak at the frequency of the island oscillations; when the island 
oscillates harmonically, the peak reduces to a J-peak. We show that knowledge of the SET transport 
properties helps to determine in what way the island oscillates, to estimate the amplitude, and the 
frequency of the oscillations. 

PACS numbers: 73.23.Hk,85.35.-p,72.70.+m 



The interplay of electric currents through nanostruc- 
tures with their mechanical degrees of freedom has at- 
tracted a lot of interest recently, both from the experi- 
mental and theoretical side [lIllllllllESIEIi- One 
of the central questions of this new field of nanophysics is 
how the vibrations of the oscillating part of a nanodevice 
influence its transport properties and vice versa. A num- 
ber of nanomechanical devices were investigated in the 
last years, e.g., so-called single-electron shuttles 
On the theoretical side, it was shown recently that elec- 
tric currents passing through a dirty nanowire can stimu- 
late its vibrations ^10..] . Indications for thermal vibrations 
of suspended single-wall nanotubes doubly clamped be- 
tween two contacts were observed 

Recently, the nanomechanical properties of single- 
electron transistors (SETs) have started to attract at- 
tention. SETs built from carbon nanotubes are a natural 
candidate for this type of question. For instance, it was 
shown that the equilibrium shape of a suspended nan- 
otube studied as a function of a gate voltage shows fea- 
tures related to single-electron electronics, e.g., Coulomb 
"quantization" of the nanotube displacement 7]. 

In this Letter, we discuss how vibrations of the central 
island of the SET change the current and the noise. We 
show that the transport characteristics of the SET differ 
for islands oscillating thermally or harmonically. The 
current noise spectrum has a peak at the frequency of 
the island oscillations that reduces to a S-peak when the 
island oscillates harmonically. Therefore, knowing the 
transport properties of the SET can help to determine in 
what way the island oscillates, and to find the amplitude 
and frequency of the oscillations. 

The system that we want to study - a SET with a 
mechanically oscillating island - is sketched in Fig. ^ 
We assume that the island is coupled to the left (L) and 
right (R) leads by tunnel junctions but can mechanically 
vibrate. A suspended nanotube is a possible experimen- 



tal realization of this scenario |7|. 

The charge of the island is coupled to the leads and 
the gate through the capacitances Cl, Cr, and Cg{z); 
z is the transverse deviation of the island from its equi- 
librium position. When the island oscillates, the gate 
capacitance changes with the position of the island cen- 
ter z, and therefore the transport properties of the SET 
change. 




FIG. 1: (a) Sketch of a single-electron transistor (SET) with 
an oscillating island. The island is coupled to the left (L) 
and right (R) leads by tunnel junctions, and its capacitance 
to the gate Cg{z) depends on the coordinate z that measures 
the deviation of the island from its equilibrium position, (b) 
Equivalent circuit of the device. 



We assume that electronic transport through the SET 
can be described by sequential tunneling. In this case, 
it is governed by four tunneling rates [ij, : the 
rate for electrons to tunnel onto the central region from 
the left (r^^jj_|_;^) and right (r^_^„_,_j) and the rates for 
electrons to tunnel off the central region. The rates can 
be calculated via Fermi's golden rule. The energy change 
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corresponding to the first tunnel process is 

e{Vg-VL)+en, 



e 



1 



CR{VH~Vg)+CLiVL-Vg[ 



where = Cl + Cr + Cg{z). Defining 7(e) = -e/(l 
exp(/3e)), the rates can be written as 



e'RL 



+n+l 



many electrons go through the SET and thus it is pos- 
sible to neglect the time derivative in the master equa- 

(1) tion (01. Then all the methods used for current calcu- 
lation in "usual" SETs are applicable to the case with 

' the oscillating island 0|. Averaging the current over 

time can be replaced by averaging over P{z), the den- 
sity of the probability distribution for the deviation z, 
i.e., I — J P{z)I{z)dz. If the island oscillations are 
thermally activated, P{z) oc exp(— z^/2(z^)t), where 

(2) = kBT/mujl. 



where Rl is the resistance of the left junction. The other 
rates can be written similarly. The charge state of the 
SET is characterized by the probability Pn{i) to find n 
excess electrons on the island. The time evolution of 
Pn [t) is governed by the master equation 0, |3 [3| : 



iPi-\-a — •n-\-aPi ) 



(3) 



a=±l 



where ri_+,- — Ffl^,, -I- Ffl^,-. The current and all its cu- 
mulants can be expressed through T and p |16l | . 

Equation Qhas to be supplemented with the equation 
describing the oscillations of the island. The electrostatic 
force that acts on the island is f{z) = dzCg{Vg — 

2/2, where ^ = {ClVl + CrVr + CgVg + 9)/Ce is 
the potential of the island, and q — ^„ Pn^e the average 
charge. Then the motion of the island can be described 
by a Langevin equation, 



z + + cj^z = y + [/(z) - /(0)]/m 



(4) 



where m is the island mass; 77 ^ oJo/Q where Q is 
the quality factor; finally, y is a random force simulat- 
ing the interaction with the thermal bath, {y{t)y(t')) = 
2rimTS{t - t') [3. Typically Q - 10^ - 10* Q and 
7? wq. If the island is a nanotube, Cg is a functional 
of the deviation z{x) where x is the coordinate along the 
nanotube The rates depend only on integral quan- 
tities like / z{x)dx/L {L is the length of the nanotube). 
Their dynamics can be described by Eq. 10} if the am- 
plitude of the nanotube oscillations does not exceed its 
diameter by several orders of magnitude. 

The current going through the left junction 0| is 



(5) 



We are interested in the current averaged over a time 
interval r much larger than the characteristic period Tq 
of the island oscillations: / = J^^ I{t, z{t))dt/2T, r ^ 00. 

The typical frequency of micromechanical oscillations 
is ujo ~ lOOMHz. If electrons tunnel through the SET 
with a similar frequency, the current will be of the or- 
der of / ^ eujo ^ 10~^^A. However, the current going 
through a typical SET is usually several orders of mag- 
nitude bigger. Hence, during the oscillation period Tq 



If the island oscillates harmoni- 
cally, z{t) = zosin(a;i), then P{z) = I/tt^^/i 
In these expressions, the driving terms 



f in Eq. ^ 

that couple the current in the SET with its mechanical 
degrees of freedom were neglected. This term is usually 
much smaller than ujqZ on the left-hand side of Eq. Q 
(e.g., for the SET parameters in recent experiments, see 
Ref.'?!): the small parameter is zq max^ dz InC^. The driv- 
ing terms may become important, e.g., when an a.c. bias 
near the resonance frequency ujq is applied to the termi- 
nals) of the SET. Then the a.c. bias drives the island 
oscillations and the /-terms in Eq. Q are important to 
determine the amplitude zq of the oscillations 

In general, the current and noise in a SET cannot be 
calculated analytically for arbitrary transport voltages 
[ifij i even when the island is static. Analytical progress 
can be made if we restrict ourself to the case of small 
driving voltages near the onset voltage, and temperatures 
much below the charging energy e^/Cs, i.e., 7(e) w ^(e) 
in Eq. In the case of a static island, the performance 
of the SET as a transistor and electrometer reaches an 
optimum in this regime ^20]. In this region, the trans- 
port characteristics of the SET are also most sensitive 
to mechanical oscillations of the island, so this regime 
is the most interesting. Only two states of the island 
have to be taken into account; the probability p has only 
two nonzero values pm Pn+i IIS- If Vl < Vr an elec- 
tron enters the island with the rate ^n^n+ii^) from the 
left lead and goes away with the rate ^n+i^ni^) into the 
right lead. The average current will be 



/ = / dzP{z)- 



n+1- 



(6) 



where n = [-CfVg/e], C^°^ ee Cg{z = 0), and [. . .] 
means the integer part. Assuming that the capacitances 
depend only weakly on z, Eq. © can be expanded with 
respect to z. To proceed, we define 



Jiz) = - 



e Rr/\E^^^_,, + Rl^ER^^^.^ 



(7) 



Using AE^'^j^^^^ = eiVR - Vl) - AE^^^_^_^ and defining 
zi, Z2 to be the roots of the equations AE^_^^_^^ — and 
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FIG. 2: Current gate-charge characteristics for symmetric 
SET {Rl = Rr) with Cg > Cl.r, Vl = -Vr = V/2, 
Qo = -CfVg. (a) V = 0.2, (|V|/2 < Ksc), (b) V = 0.5, 
{\V\/2 « Ksc), (c) 1/ = 1, (|V|/2 > Ksc). The voltage is 
measured in units of \e\/2Cs{z = 0), the current in units of 
Jo = (e/Cf>)/(i?L + Rr). The dashed curve corresponds 
to a static island, the solid curve to a harmonically oscil- 
lating island, z — zosm{uJot); zo{dzCg)/Cg — 5.6 ■ 10^^ 
(this is typical for SETs where the island is a nanotube 
[3; then zo ~ 5r, where r is a typical nanotube diame- 
ter). The dotted curve in Fig. ^ illustrates what happens 
if V;sc/(|e|/2C^"') = 5 > 1 and Eq. © is not valid. The 
dash-dotted curves correspond to the case of thermal mo- 
tion; the thermal average {z^)t = kBT/muio is chosen to 
be equal {{zosm{ujot))^)t — Zq/2. The integer part of Qo/e 
is the number of electrons on the island in the static regime 
when Vl = Vr = 0. The curves are periodic in the static case, 
but not if the island oscillates. The areas under the peaks in 
the static and in the dynamical cases are the same. 



Ai;^+i^„ = 0, we get from Eq. ® 



^max(zi,Z2) 

/« / dzPiz)[JiO) + 

J min(^i ,Z2) 

z{d.J{z)U^o) + z'l{dU{z)U^o)] ■ (8) 



This formula is valid also for Vl > V/?. Using Eq. 
we find Z2 - « eiVR - Vl) / dzen\z^o- Thus, if zq is 
a characteristic amplitude of island oscillations then it is 
natural to define the voltage scale 



Vosc — —{dz£n\z^o) 

e 



(9) 
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FIG. 3: Coulomb staircase, Rr/Rl = 1000, Vl = -Vr = 



V/2, Cl = Cr = Cf \ Qo = 40.6 {C^"' and other parameters 
are as in Fig. |5J . The average square amplitude of the island 
deviations is equal to that of the harmonically oscillating is- 
land. The Coulomb steps are smeared on the scale Vosc- For 
other choices of Vg the staircase graph looks similarly. 



(0) 



In the hmiting case Cl , Cr <C C, 



(0) 



zodzCg 



-e{n + l/2) 



(10) 



Eq. ||SJ) is valid if V^sc < e/Cs- If the driving voltages 
applied to the SET terminals are much larger than Vosc, 
the integration limits in Eq. can be extended to in- 
finity because they far exceed zq, the scale of decay of 
P{z). The second term in Eq. © vanishes and 



^" = ^(--0) + i(z^)^/(z)|. 



(11) 



where (z^) = / P{z)z'^dz. The first term in Eq. ^fTT}i is 
the current for a static island. If the driving voltages 
applied to SET terminals are smaller than Vq^c then the 
second term (linear in z) in Eq. (jSJ does not vanish; in 
this regime the current-voltage characteristics is largely 
influenced by island oscillations. The small parameter 
in the expansion Eq. ((HJ is zodzCg. The second term 
in Eq. Hll|l is of second order in this parameter. If we 
investigate the I — V characteristics, island oscillations 
will smear the Coulomb gap within a voltage band of 
width of order Vosc- 

The I — Vg characteristics of a symmetric SET [Rl — 
Rr, Cl = Cr <^ Cg) with an oscillating island is shown 
in Fig. 13 The dashed curves correspond to the case of 
the static island. The solid and dash-dotted curves show 
the case of a harmonically oscillating island or an island 
subject to thermal equilibrium fluctuations. Figures [^a- 
c) also show how the I — Vg characteristics change when 
the driving voltage V is smaller, of the order of, or larger 
than Vosc- Within each peak of the curves, n is constant. 
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FIG. 4: Fano factor for V < Vosc The dashed hne corre- 
sponds to the static case, the sohd hne to the harmonicaUy 
osciUating island, and the dash-dotted line to the thermally 
driven island. The average square amplitude of the island os- 
cillations are equal for the harmonically oscillating island and 
the case of thermal equilibrium. 

therefore V^sc is also constant within the peak. The most 
interesting case is shown in the first panel of Fig. |2Ia). 
For the harmonically oscillating island, the peaks split 
and their width becomes larger with the characteristic 
scale Vosc', when the island moves due to thermal activa- 
tion the peaks do not split but their width also becomes 
larger with Vosc ■ Thus, the type of motion of the island 
leaves a characteristic trace in the I — Vg plot. Equa- 
tion © well describes I — Vg characteristics when the 
peaks do not overlap, like in Figs. and b. It follows 
from Eq. lO that the areas under the peaks the in static 
and dynamical cases are equal. The I — Vg characteristics 
is periodic in the static case, but not periodic for an os- 
cillating island because Vi^sc changes from peak to peak, 
see, e.g., Eq. ((Tn|) . For this graph, we used the parame- 
ters of the nanotubc model of the island (see Ref. 0); we 
chose, e.g, Cg{z) — L/2ln{2{R—z)/r), where r = 0.65nm 
and L = 500nm are the nanotube radius and length, and 
R = lOOnm is the distance to the gate. If Cl, Ch and Cg 
are of the same order, the I—Vg characteristics is qualita- 
tively similar to what is shown in Fig. [5] If the junction 
is asymmetric, Rr ^ i?L, the peaks in Fig. [21 become 
nonsymmetric. Figure 13 shows the I — V characteristics 
for Rn ^ Rl. The oscillations of the island smear the 
Coulomb gap within a voltage range of the order of Ksc- 
We will now discuss the current noise in a SET with a 
moving island. When the island oscillates, the irreducible 
current-current correlator 5 (r, 0) = ((/l(6+t/2)/l(0 — 
r/2))) depends on both r and Q (rather than only on r 
as in the case of the static island). However, since the 
charging events in the SET are correlated on time scales 
much shorter than the period of the island oscillations, 
the dependence of 5 on t is much stronger than on Q, 
and the zero-frequency noise can be found as / S{t, Q)dT, 
where the bar means averaging over 0. In other words, 
the low-frequency noise can be calculated at a given po- 



sition of the island (e.g., as it is done in Ref. H^ and then 
averaging over time as it was already done for the current 
above. The result of this procedure is presented in Fig.^ 
which shows the dependence of the Fano factor [ij on 
Vg. Here we assumed that the driving voltage is smaller 
than Vosc, i-e., the system is in the region in which the 
influence of the oscillations of the island on the trans- 
port properties of the SET is maximal. The Fano factor 
dips get split due to harmonic oscillations of the island; 
the scale of the splitting is Vosc- In contrast to that, the 
dips are washed out by thermal equilibrium oscillations 
of the island. Here, we assumed that the noise frequency 
is much below the frequency of the island oscillations, luq. 

We now consider the noise spectral density, i.e., the 
Fourier transform S{uj,Q) of S{t,Q). If w approaches 
luq, the correction to the noise from the motion of the is- 
land is ^ (dzl)'^ {{z'^))uj, where {{z'^))uj is the irreducible 
correlator of the deviation z at frequency to. It has a 
(5-peak at ojq if the island oscillates harmonically. In 
contrast, if the island moves due to thermal activation, 
{{z^))^ = 2r]kBT/m{{uj^ - wg)^ + to^r]^) M, the noise 
peak has a width of the order of the oscillation damp- 
ing factor 77, sec Eq. A similar result was obtained 
numerically in Ref. 22. Thus, measuring the noise spec- 
trum allows to find the frequency of island oscillations 
and gives information on the nature of the oscillations. 

In conclusion, we have discussed how vibrations of the 
island in a SET change its transport properties. The 
transport characteristics of the SET can be used to de- 
termine the nature of island motion, in particular, to 
estimate the amplitude and frequency of its oscillations. 
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